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We study the generation of primordial fluctuations in pure de Sitter inflation where the quantum 
scalar field dynamics are governed by polymer (not Schrodinger) quantization. This quantization 
scheme is related to, but distinct from, the structures employed in Loop Quantum Gravity; and it 
modifies Standard results above a polymer energy scale M*. We recover the scale invariant Harrison 
Zel'dovich spectrum for modes that have wavelengths bigger than M^ 1 at the start of inflation. 
The primordial spectrum for modes with initial wavelengths smaller than M* 1 exhibits oscillations 
superimposed on the Standard result. The amplitude of these oscillations is proportional to the 
ratio of the inflationary Hubble parameter H to the polymer energy scale. For reasonable choices 
of M*, we find that polymer effects are likely unobservable in CMB angular power spectra due to 
cosmic variance uncertainty, but future probes of baryon acoustic oscillations may be able to directly 
constrain the ratio H/M+. 



I. INTRODUCTION 

It is commonly believed that quantum gravity ef- 
fects may significantly alter "standard" physics near the 
Planck scale. For example, string theory posits that ex- 
tra dimensions with compact topology will become vis- 
ible at energies approaching Mpi, while loop quantum 
gravity asserts that continuous classical spacetime is re- 
placed by quantum spin networks on small scales. Unfor- 
tunately, the huge discrepancy between the Planck scale 
and typical energies in the nearby universe make it virtu- 
ally impossible to experimentally or observationally test 
such ideas. In fact, the only available data comes in- 
directly through measurements of the cosmic microwave 
background: If one accepts the inflationary paradigm, 
this thermal relic of the hot big bang depends on the 
spectrum of primordial perturbations generated when the 
temperature of the universe was a few orders of magni- 
tude less than the Planck scale. This makes inflation the 
only known phenomenon which both involves Planckian 
energies and has measurable consequences for the observ- 
able universe. It is therefore a crucial issue for theories 
of quantum gravity to predict how "new physics" near 
the Planck scale affects the spectrum of primordial per- 
turbations generated from inflation. 

A useful way of thinking about the effects of new 
physics in inflation involves viewing the quantum gen- 
eration of inflationary fluctuations as a trans-Planckian 
problem [T]. The idea is as follows: if one takes the 
very largest scale cosmological perturbations which are 
relevant for observations and tracks them backwards in 
time, one finds that at some finite time during the in- 
flationary epoch their physical wavelengths will become 
smaller than the Planck length. Hence, the early time 
evolution of such modes will necessarily be sensitive to 
any new physics manifest at small scales, which then im- 
plics that thcrc should be some imprint of very high en- 
ergy phenomena on very large cosmological distances. Of 
course, the key questions are the amplitude and nature 



of these effects, which in turn depend on the nature of 
the small scale modification. 

In the literature, various authors have considered ad 
hoc modifications to scalar field dispersion relations [I]- 
0], due to non-commutativity [5HE], or modified uncer- 
tainty relations [5HT3"] . There have been attempts also to 
calculate trans-Planckian contributions to the primordial 
power spectrum in a model-independent way by impos- 
ing initial conditions on a "new physics hyeprsurfacc" 
[T3HH1- Recently, effects arising in Horava-Lifshitz grav- 
ity have been reported [18]. A feature of many (but 
not all) of these studies is that short distance effects su- 
perimpose oscillations on the conventional scale- invariant 
power spectrum with amplitude (H/M+) 1 , where M* is 
the energy threshold above which the modifications are 
important, H is the inflationary Hubble parameter, and 
the power 7 > 1 depends on the model. 

In this paper, we explore a diffcrent class of "new 
physics" suggested by the "background independent" (or 
"polymer" ) approach to quantization that is deployed in 
loop quantum gravity (LQG) |19j . In this programme, 
classical geometric variables such as the metric are rep- 
resented at the quantum level by graphs on spatial 3- 
manifolds known as spin networks. Fundamental geo- 
metric information such as areas, volumes and their evo- 
lution are encoded by densitized triads and holonomies of 
conncction 1-forms over the edges of these graphs. The 
key point is that while quantum operators corresponding 
to the holonomies are well-defined, operators correspond- 
ing to the connection 1-forms themselves are not. This 
implies the Hilbert space of LQG has distinct properties 
from the Standard one underlying Schrodinger quantum 
mcchanics and quantum field theory. 

The novel features of this approach to quantization 
are best illustrated by considering ordinary quantum me- 
chanics [20] : Consider a particle moving in one dimcn- 
sion and described by a position x and its conjugate mo- 
mentum p. In the conventional Schrodinger quantization 
(SQ) of the system the partiele's state is described by 
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an element of a Hilbert space in which the action of the 
position x and momentum operators p are well-defined. 
On the other hand, an LQG-inspired quantization makes 
use of an alternative Hilbert space where the operator 
corresponding to p is not dcfincd, but the associated 
"holonomy" operator is. Since p is the generator of in- 
finitesimal translations, the appropriate identihcation of 
its holonomy is the finite translation operator U\ (whose 
action is to displace the particle by a distance Ax = — A). 
The quantization algorithm based on this Hilbert space is 
callcd "polymer quantization" (PQ) since it is motivated 
by the spin network structure of LQG, where the excita- 
tions of the gravitational ficld occur along the edges of a 
graph; i. e. they are one-dimensional like a polymer. 1 

The lack of a natural momentum operator in PQ may 
seem alarming, but one can easily define an effective p 
by constructing a simple finite diffcrcncc stencil for the 
Schrodinger momentum operator id x using finite trans- 
lations U\ it . The charaeteristie size of this stencil A* is 
an arbitrary fixed parameter of the quantization that de- 
fines a polymer energy scale M*. We expect to recover 
ordinary SQ at energies less than M+ (since our finite 
difference approximation to p will be very good in that 
regime) while at higher energies we would expect the pre- 
dietions of PQ and SQ to diffcr substantially. 

This expectation has been explicitly confirmed by cal- 
culating the spectrum of a polymer quantized simple har- 
monic oscillator of mass m [20] [23] 121] • One finds that 
the energies E n of eigenstates of the Hamiltonian ap- 
proximate the well-known SQ values when mE n /M+ <C 
1. Similar results are available for the Coulomb [25], 
inverse-square [25] , and other [27] [25] spherical poten- 
tials; though the issue of boundary conditions at the ori- 
gin of spherical coordinates must be handled carcfully 
[25] . The limit in which one can obtain reasonable poly- 
mer approximations to Schrodinger wavefunctions has 
also been considered [30] . 

Polymer effects have also been studied extensively in 
the context of quantum cosmology 31-33J. Specifically, 
the polymer treatment of geometric quantities in FRW 
models preserves the predietions of general relativity at 
low curvature while replacing the big bang singularity 
with a big bounce when the density of the universe is 
~ [34] . Conversely, polymer quantization of a ho- 
mogeneous and massless scalar in the early universe has 
been shown to replace the big bang with a past eternal de 
Sitter phase with Hubble parameter H ~ M* /Mp\ [35] . 

It is natural to try to extend these results from poly- 
mer quantum mechanies to quantum field theory. In that 



vein, the PQ of a scalar field in Minkowski space has been 
considered [36], assuming compaet topology [37], using 
semi-classical approximations |38j , via an effective spatial 
lattice |39j , and by direct quantization of Fourier modes 
|24) . In the last approach, it was shown that Fourier 
modes with e ik ' x spatial dependence exhibit exotic poly- 
mer behaviour if |k| 3> M*. 

In this paper, we use the techniques introduced in |24j 
to study the PQ of a scalar field in a de Sitter inflationary 
universe. The motivation is obvious: since the physical 
wavenumber of a given Fourier mode is inversely propor- 
tional to the scale faetor in an expanding universe, its 
behaviour will be dominated by polymer effects in the 
asymptotic past. Hence, we would expect that the PQ 
of a scalar field during inflation will result in potentially 
observable modifications to the primordial perturbation 
spectrum. The current work confirms and quantifies this 
expectation. 

The organization of the paper is as follows: In j[TT] we 
recall the Standard textbook caleulation of the primor- 
dial power spectrum as well as an alternative formula- 
tion based on quantization of individual Fourier modes. 



In [III we deseribe how mode-by-mode quantization is 
achieved in the Standard Schrodinger picture, while in 
§IV| we present the caleulation in the polymer formal- 
ism. In SjV] we present numerical and semi-analytic re- 
sults for the polymer primordial spectrum, and in §VI| 
compare them to observations of the cosmic microwave 
background and large scale structure of the universe. 
We summarize and diseuss our main results in f VII 



1 It is important to note that PQ has been proposed as an alter- 
native to Standard quantum theory independently of any quan- 
tum gravity considerations: In particular, Halvorson 12 1 1 pro- 
posed such a quantization as an alternative to SQ that allowcd 
for normalizable position eigenstates in the Hilbert space. The 
existence of such states implies that the Standard uncertainty 
principle does not hold in the polymer picture 1221 . 



The appendices give a technical introduetion to polymer 
quantum mechanies j |A"] list cosmological scaling relations 
uscd throughout the paper SjB] and derive some technical 
formulac SjCl 



II. GENERATION OF PRIMORDIAL 
PERTURBATIONS IN A DE SITTER UNIVERSE 



In this seetion, we review the caleulation of the spec- 
trum of primordial perturbations in a de Sitter inflation- 
ary universe using two complementary methods: The 
flrst is based on the quantization of the scalar field in 
real space and the subsequent Fourier decomposition of 
the quantum operators. The generation of fluetuations 
follows from the faet that quantum operators obey the 
classical equations of motion; i. e.; this approach uses the 
Heisenberg picture. The second method involves flrst 
Fourier decomposing the field, which reduces the system 
to a collection of independent oscillators with time de- 
pendent parameters, and then quantizing each oscillator. 
In this case, the generation of fluetuations follows from 
the solution of the resulting one-dimension Schrodinger 
cquation with a time dependent mass and potential for 
the wavefunction (i. e., this approach uses the Schrodinger 
picture). We will make use of the latter approach when 
considering the polymer quantization of the scalar field. 
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A. Quantization in real space 

We consider a massless scalar field propagating in a de 
Sitter background: 2 

ds 2 = -dt 2 + a 2 (i)dx 2 , a(t) = expHt. (1) 

The Hamiltonian of the scalar field is given by 



Mj, = I d 3 xa 3 



1 

2^ 



1 



{v<t>y 



(2) 



Here, tt is the momentum conjugate to 4> such that 
{</>(£, x), 7r(£, y)} = (5^ 3 ^(x — y). To quantize this system, 
one customarily promotes both <f> and 7r to operators and 
imposes the commutation relations [</>(£, x), ir(t, y)] = 
«<!>( 3 )(x — y). Then, the field operator is decomposed into 
Fourier modes via 



- ikx 



ak + / k (*) e 



+ik-x-t 
"k 



(3) 



Minkowski vacuum state is recovered in the appropriate 
limit of the parameter space, which picks out the solution 



./k 



H 

TIP 



Ha 



D ik j Ha 



(7) 



This choice is referred to as the Bunch-Davies or adia- 
batic vacuum, and it gives rise to the familiar scale invari- 
ant Harrison-Zel'dovich (HZ) spectrum Vuz = (H/2tt) 2 . 



B. Quantization in Fourier space 

The algorithm we just described involved quantization 
of <j) first and then decomposition into Fourier modes. 
However, there is a equivalent procedure that involves 
Fourier decomposition and then quantization. The first 
step is writing 



<K*,x) = 4^I>k(i)< 

V ^0 ,. 



d A x e- lk ' x 0(i,x) 



(8) 



where Vq is the fiducial volume (with respect to the flat 
3-metric diag{\, f, f)) used in our box normalization i. e. 



Vo 



d A x. 



(4) 



In the Heisenberg picture, <j> satisfies the classical equa- 
tion of motion U(f> = 0, which implies that 



/k + 3tf/k+^/ k = 0. 



(5) 



The ladder operators satisfy the commutation relation 
[etki^k] = 1- During inflation, we assume that the field 
is in the vacuum state annihilated by the Sk operators; 
i. e. &k IV 1 ) = 0- Then, the power spectrum of primordial 
perturbations generated during inflation is 



mu) 



2tt 2 ' 



<0k> = WklV>) = l/k| 2 . 



fc<aff 



(6) 

Hence, we see that the power spectrum is entirely deter- 
mined by solution of the classical equation of motion 
Now, there are infinitely many solutions of the ODE and 
each particular solution will pick out a different quantum 
vacuum state; i. e., there is an ambiguity in the calcula- 
tion. This is commonly resolved by demanding that the 



2 Recall that in a de Sitter background, the equation of state 
p + p = implies that metric perturbations are decoupled from 
scalar field fluctuations, so (f> is automatically a gauge-invariant 
quantity. 



with a similar expansion for Tr(t, x). After a suitable 
rcdefinition of the indcpendent modes to ensure that <fi 
and the rcdcfincd </>k's are real, the Hamiltonian is 



^ = 22^ = 22 



2a 3 



2a 



(9) 



with the Poisson bracket {<f>k,TTk>} = <5k,k'- 

This expression of the Hamiltonian ^ implies that 
we can view the classical system as a collection of inde- 
pendent oscillators with time-dependent parameters and 
labeled by the wavevector k. The quantum state of the 
field is of the form 



'l^k 



(10) 



where each of the l^k) satisfy the time-dependent 
Schrodinger equation (TDSE), 



J^klV'k) = id t \i^)< 



(H) 



As described by [13 ST] , to determine the spectrum of 
primordial perturbations generated during inflation one 
first solves the (TDSE) for the "ground state" |^ k ) = |0 k ) 
of each mode. Then, the power spectrum is given by: 



= (0k|^|0 k ). (12) 



fc<_ffa 



This mcthod is useful because it reduces the full quantum 
field theory to problem in quantum mechanics. This is 
crucial, because PQ is much easier to deal with a purely 
quantum mechanical setting. 
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III. SCHRODINGER QUANTIZATION OF AN 
INDIVIDUAL FOURIER MODE 

We are ultimately interested in the solution of the 



TDSE ( 11 ) with the polymer representation of the Hamil- 
tonian operator governing a single Fourier mode, but as a 
prelude we review how the calculation works using Stan- 
dard Schrodinger quantization (SQ). Rccall that we can 
represent arbitrary quantum states as wavefunctions de- 
pending on either "position" or "momentum" , which in 
our case correspond to </>k or 7Tk, respectively. In the case 
of SQ, both choices are very similar due to the symmetric 
form of the simple harmonic oscillator Hamiltonian 
However as described in Appendix[Aj the polymer quan- 
tization of position and momentum operator are handled 
in quite different ways, which has the net effect of making 
it easier to work with momentum space wavefunctions. 
Hence, we assume arbitrary quantum states are repre- 
sented by functions of tt^: 



(TTklV'k) = 1p(t,TTk), 



(13) 



where |7Tk) is a momentum eigenstate. The action of cf>^ 
and 7Tk on these wavefunctions is simply: 



(TTkl^klV'k) 
(TTk^kl^k) 



(14a) 
(14b) 



It is easy to confirm that this operator representation 
respects the commutation relation [</>k,7Tk] = i- Using 
these, we find that the action of the Hamiltonian operator 
on an arbitrary state is: 



(7T k |<«k> 



1 

— 'i 

2/z 



\iuj 2 d 2 



^(t,7Tk); (15) 



where we have defined the time-dependant parameters 



/t : 



= k/a. 



(16) 



With the representation (|15j), the time-dependant 
Schrodinger equation, 



(17) 



reduces to a PDE for the wavefunction ip. If we now 
change to the conformal time r\ and a dimensionless mo- 
mentum y, 



1 

Ha' 



y 



7Tk 

JJIuj 



= —krj 




and re-scale the wavefunction as 



we find that 



V- k v^(v,y) 



exp 



2kriJ 



2<V 



(18) 

> (19) 

(20a) 
(20b) 



We see that the new wavefunction is just a solution of 
the familiar TDSE for an oscillator with unit mass and 
frequency — and the new effective Hamilt onia n H is time 
independent. 3 Normalizable solutions of (20) satisfy the 
boundary conditions 



*(77,±oo) = 0, 

and are well known: 

oo 

*( ?7 ,jy) = ^ C „* Il (y)e-^ +1 / 2 ) fe ", 
with 



(21) 



n=0 



n=0 



* B (») = ( x /7r2"n!)- 1 / 2 J ff„( l/ ) e -^/ 2 . 



(22) 



(23) 



Here, H n is the Hermite polynomial of order n. Notice 
that the constants in the 'J/ definition ensure that if \P is 
normalized with respect to integration over y, ip is also 
normalized with respect to integration over w^: 



1 = 



dy\^(f],y)\ 2 = / dir k \i>(t,Tr k )\ 2 



(24) 



Having now obtained the explicit solution for the 
TDSE, we need to identify the ground state wavefunc- 
tion. The most natural choice is to use the ground state 
of H, 



tffoy) = ^ (y)e- lfc "/2 = n -l/i e -{y 2 +ikv)/2 } 



(25) 



since this state minimizes the expectation value of the 
effective Hamiltonian for all r]. Indeed, it is precisely this 
choice that reproduces the familiar Bunch-Davies result: 



4) 





ip* ( 


J oo 




H 2 




k 3 . 


f — OO 


H 2 




2P 





. d 

y 2 \*\ 2 

k 2 
H 2 a 2 



Ip d7Tk 
k 2 



H 2 a 2 1 y 



dy 



(26) 



Just as in the Standard calculation, we would have ob- 
tained a different answer had we imposed different con- 
ditions on the quantum state of the system. 



3 We stress that the effective Hamiltonian H should be viewed as 
a linear differential operator rather than a true quantum Hamil- 
tonian; for example, it does not generate the time evolution of 
observables. Having said this, we will use the terms "ground 
state", "energy eigenstate", "eigenenergy" , etc. when referring 
to solutions of the ODE eigenvalue problem for H in this subsec- 
tion and below. 
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IV. POLYMER QUANTIZATION OF AN 
INDIVIDUAL FOURIER MODE 

A. Formal solution of time dependent Schrodinger 
equation 

We now seek to find a "ground state solution" of the 
polymer version of the TDSE governing the amplitude of 
a given mode in k-space (11), from which we can calcu- 
late V (f,- We again work with momentum wavefunctions 
ip(t, 7Tk). Our treatment will largely follow the discus- 
sion of the polymer quantization of a particle moving in 
a one-dimensional potential in Appendix [XJ with a few 
notable exceptions. 

As in the Schrodinger case of we represent the 
quantum state of an individual Fourier mode in a basis 
kk>: 



(TklV'k) = VK^k) 



(27) 



However in the polymer scenario, |7Tk) is not interpreted 
as a momentum eigenstate because the momentum op- 
erator does not exist. But we can define an operator 
U\ that will be seen to correspond to finite translations 
of the held amplitude: 



U\\~>Pk) = exp(£Aa 3/2 7r k )^(t,7r k ). 



(28) 



This definition is similar to the one given in Appendix[A} 
except for the a~ 3 / 2 faetor in the argument of the expo- 
nential. As in Ref. [35 , this has been included to ensure 
that U\ transforms as a scalar under the dilation x — > £x 
(c/. Appendix |b|) , which will ensure we recover the cor- 
rect scaling of the effective momentum given below. We 
define the operator corresponding to the field amplitude 
4>k in the same manner as the SQ case: 



(29) 



It is easy to confirm the following commutation relation 
holds: 



[<Pk,U\] = -AL> A 



(30) 



where we have defined the smeared Fourier amplitude 
operator 



<Pk = a 3/2 ^k, 



(31) 



which transforms as a scalar under x — > £x. This allows 
us to further interpret U\: Suppose |</? k ) is an eigenstate 

of Ifik- 



(32) 



Then Ux\^k) will be an eigenstate of <£>k with eigenvalue 
tpk - A: 



Vk 



= (<^k - A) (U\\<Pk) 



(33) 



Hence, U\ has the effect of inducing translations of mag- 
nitude A in the smeared Fourier amplitude <^>k- 

Once we have a representation of U\, we realize the 
momentum contained in the Hamiltonian operator 
as i\ finite difference operator: 



TTk H» 7T k = 



where A* = M* is a fixed parameter with dimensions 
of (mass) -1 / 2 . Since U\ t transforms as a scalar, we see 
that 7r£ transforms like ir^ under dilations. Furthermore, 
we recover the SQ momentum operator in the appropri- 
ate limit: 

lim (7r k |7!k|V'k) = lim (^kl^kl^k) = 7TkV>(*> 7r k ). (35) 

A*— ^0 a— >-oo 

It is fairly easy to confirm that this gives the aetion of 
J^k on an arbitrary state as 



a 3 / 2 
2i\± 



(34) 



(TTkl^klV'k) 



sin 2 (A7r k ) d 2 
2/xA 2 



2 dnl 



^(t,7T k ), (36) 



where /i and u are defined as above (16), while 
A = Aa~ 3 / 2 . 



(37) 



In the A-7Tk — > limit, we see that (36) reduces to the 



SQ expression (15). As deseribed in detail in Appendbc 



|A| if we restrict ourselves to one super-selected seetor of 
the polymer Hilbert space the appropriate inner produet 
between states is: 



7T/2A 



£*(t,ick)ip(t, 7r k )(i7rk 



(38) 



'-7T/2A 

Furthermore, wavefunctions will satisfy 



ip(t, -7T/2A) = exp(in^./A)ip(t, tt/2A), (39) 

where 0k/A G [0, 2) is a constant (this is the lattice offset 
of the super-selected seetor, as diseussed in Appendix[A|. 

As in the SQ case, our task is to now solve the TDSE 
with the polymer Hamiltonian. We again change coordi- 
nates and re-scale the wavefunction as in Eqs. ( 18 1 and 



( 19 1 to obtain 



d 



1 d 2 
2<V 



H = k 



s\VL 2 {g 1 /' 2 y) 

2g 



(40a) 
(40b) 



4 It is important to note that this realization of the momentum 
operator in polymer quantum mechanies is not unique; i. e., it 
represents a new type of quantization ambiguity in addition to 
ones already present in conventional Schrodinger quantization. 
Stated another way: our quantization seheme is defined by our 
choice of fundamental operators ip^ and U\ as well as our spec- 
ification of the momentum operator ( | 34 [ > . Whether or not other 
finite difference representations of the momentum yield different 
physical results is an open qucstion. 



G 



where the polymer coupling is defined by 



(41) 



Unlikc the SQ case, the efFective Hamiltonian H is a func- 
tion of time via its dependence on the coupling g. At late 
times rj — > and g — > 0, which gives H — > \{y 2 — d 2 ,); i. e. 
we recover the SQ Hamiltonian. Note that under this 
change of coordinates and scaling of the wavefunction, 
the inner product becomes 

(6c|V<k) = j Z*(r 1 ,y)*( V ,y)dy, (42) 

yel 

where ^(t, 7Tk) and £(i, 7Tk) are the images of ^/(j?, y) and 
5(7/ , y) under the wavefunction transformation defined by 
(19), respectively, and 



We now seek a solution of (40) in terms of the eigen- 



functions of the time-dependent Hamiltonian H. To fix 
these eigenfunctions, we need to specify boundary con- 
ditions on ty(r),y). The relationship (39) and the wave- 
function transformation ( 19 ) imply 5 



|*(ry,-7r/2V9)| 2 = |^,7r/2^)| 5 



(44) 



Now, we impose that the evolution of an arbitrary state 
|?/>k) be unitary; that is, 



di) 



(V'klV'k 



0. 



(45) 



Carrying out the differentiation by making use of (40), 
d42l, and d44l we obtain: 



*(77,±7r/2Vff)=0; 



(46) 



that is, $ must satisfy Dirichlet boundary conditions. 
Note that this is consistent with the boundary condition 
(A20), which enforces that the spectrum of the Hamil- 



tonian in polymer quantum mcchanics is independent of 
lattice offset. Finally, note that (46) recovers the SQ 



boundary conditions (21) in the g — > limit 



The solution of the energy eigenvalue problem 

H(#„(l, y) = e„ (77)*« (»7, y), (47) 



subject to the boundary conditions ( 46 ) is given explicitly 



in terms of Mathieu (elliptic sine) functions 
^2 



*n(ff,tf) 



cj 1/4 se„ +1 



1 „ 7T 



(48) 



This could have also been deduced by demanding that the prob- 
ability amplitude l^l 2 s har e the same periodicity as the effective 
potential appearing in \40\ . Or, by demanding that 0k be self- 
adjoint under the inner product (|42|). 



with eigenvalues given by the Mathieu characteristic 
value functions: 



2gB nA 



1 



(49) 



for n — 0, 1, 2 . . . These form an instantaneous orthonor- 



mal basis for arbitrary functions on / satisfying ( 46 1 



K^m dy 



(50) 



yel 



Also, eigenfunctions with even n have even parity 
and those with odd n have odd parity: ^ n (— y) — 
(—l) ni ff n (y). We can use the asymptotic expansions of 
the Mathieu functions [42 to deduce 



1 



( 43 ) ?/) 



l~2 



H n (y)e- y2/2 
(n + 1 



gy 



9 < G n 

9 > G n 
(51) 

where G n = l/(n + l/2); i. e., we recover the Schrodinger 
energy eigenfunctions for small polymer coupling and 
simple trigonometric function for large coupling. Asymp- 
totic expansions of the Mathieu characteristic value 
function B n yield the following approximations for the 
eigenenergies: 



\\{n + l) 2 g, 3 >G r; 



(52) 



Plots of the energy eigenfunctions and eigenvalues are 
presented in figures [T] and [2j respectively. 

Since our effective Hamiltonian is time dependent, ar- 
bitrary solutions of pO] ) can be constructed from Hamil- 
tonian eigenstates if we allow the expansion coefficients 
to depend on ?y: 



Mi?) i 



where 9 n (v) satisfies 



,(»70) - 0, 



(53) 



(54) 



Here 770 is some initial time, and we use an overdot to 
indicate derivatives with respect to 77. Some straight- 
forward algebra reveals that the expansion coefficients 
satisfy an autonomous set of linear ODEs: 



-c n (n\n) 



E' 



(n\H\m) 



,i(6W-M 



(55) 



The normalization condition (n\n) = 1 in addition to the 
fact that the energy eigenfunctions are real imply that 
(n|n) = 0. Since the polymer coupling g is proportional 
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FIG. 1: Instantaneous energy eigenfunctions of the effective 
polymer Hamiltonian. At moderately small polymer coupling 
(left), these resemble the energy eigenstates of the simple har- 
monic oscillator, while at larger coupling (right) they reduce 
to trigonometric functions. 
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FIG. 2: Energy eigenvalues as a function of g of the effective 
polymer Hamiltonian 



to the conformal time r\ (41 ), we can recast this as a first 
order matrix ODE: 



dg 





CO 




«00 a oi • • • 




Ac, c = 


Cl 


, A = 


a 10 a ll 


, (56) 



exp 



yel 



.M* 

i 

H 

d 
' dg 



90 



-)d~g 



sm 2 (g 1 / 2 y) 
2ff 



^ndy, (57) 



for n ^ m and a nn — 0. Here, go is that value of the 
polymer coupling at time 770. Note that the matrix A is 
anti-Hermitian (A^ = —A) so the norm of c is conserved: 



dg 



(c f c) = 



(58) 



Also note that will be non-zero only if ^„ and 
^ m have the same parity; i. e., dj i+2fe+i = for k — 
0,±1,±2,... 



Initial conditions and the final spectrum of 
fluctuations 



A unique solution of the polymer TDSE will be char- 
acterized by the specification of initial conditions for the 
expansion coefficients at some time 770. The question is: 
what choice of initial conditions could reasonably be as- 
sociated with the vacuum state of a given Fourier mode? 
In the case of SQ, the answer was straightforward be- 
cause the effective Hamiltonian H was independent of 77, 
implying that if we prepared a given mode in its ground 
state at some initial time, it would stay in its ground 
state indefmitely. That is, the quantum evolution was 
perfectly adiabatic. This is not true in the PQ case: (55) 



tells us that there is non-trivial mode mixing. That is, if 
we prepare the system in the ground state at some early 
time, it will not be in the ground state at the end of 
inflation. 

The situation is akin to the trans-Planckian problem of 
inflationary cosmology considered by Martin and Bran- 
denberger [T]. In that work, the classical equations of 
motion of the infiaton werc modified on small scales in a 
attempt to account for quantum gravity effects. The net 
result was that one could not unambiguously identify a 
vacuum state in the asymptotic past due to particle cre- 
ation induced by the modified wave equation, which is 
precisely analogous to the mode mixing induced by poly- 
mer effects at early times in the current scenario. Hence, 
we will adopt the same prescription for initial conditions 
as employed in [T]: We assume that the field is in the 
instantaneous ground state at the beginning of inflation. 
That is, 



'(flo) = [l •••] t , .go = .9(7/0), 7/0 



1 



Ha 



(59) 



Here, ao is the scale factor at the beginning of inflation. 
Notice that the evolution equation for the coefficients 
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FIG. 3: Scale factor evolution of the physical wavelength of 
modes with large or small go relative to the polymer M^ 1 and 
Hubble H^ 1 scales 



( 56 1 is in terms of polymer coupling instcad of the con- 
formal time, so it useful to know that the initial value of 
g for a given mode 



fe 

K' 



fe* = M*a 



(60) 



Here, the pivot scale fc* is the wavenumber of a mode 
that has physical wavelength M~ Y at the beginning of 
inflation. As derived in Appendix[C] the numeric value 
of fc* is 



6 x 10- 
Mpc 



M* 

~TT 



10 16 GeV j 



-,65 



100 

1T 



1/12 



(61) 

where E inl — pV f 4 is the energy scale of inflation, N = 
ln(a en d/ao) is the number of e-folds of inflation, and G is 
the effective number of relativistic species at the end of 
inflation. 

Having now specified the quantum state at beginning 
of inflation, we note that observable quantities are di- 
rectly derived from the quantum state of the end of in- 
flation. We will restrict our discussion to modes with 
physical wavelengths much less than M~ x at the end of 
inflation; i. e, modes with 



.9cnd = e N g Q < 1. 



(62) 



Hence, the final quantum state of the system will be given 
by 



„end 



c(0). 



(63) 



We are also interested in modes that are well into the 
superhorizon regime at the end of inflation; i. e., modes 
with 



k 



Ha 



« 1. 



(64) 



Note that assuming H < M* implies that all superhori- 
zon modes at the end of inflation will automatically have 



.9ond 'C 1. Using the definition (19), we see that power 
spectrum of such modes is 



2tt 2 



(K) « H (*l 2 y 2 l*) 

(65) 

where we have assumed k/ Ha -C 1 to obtain the last 
expression. The state vector in this expression is given 
by the rj -> limit of ((531: 



fc<ga£T 



E 



end i \ 



(66) 



where |ti)sq are the energy eigenstates of the ordinary 
Schrodinger quantized simplc harmonic oscillator. 

To summarize, in order to calculate the power associ- 
ated with a mode with a given value of go = fc/fc* at the 
end of inflation, we must solve the matrix equation (|56|) 



subject to the initial condition (59) for the final values 
of the expansion coefficients c (0)7^ Thcn, these expan- 
sion coefficients can be inserted into the formula ( 65 ) to 
obtain V$(k). 



C. Perturbative solutions 

1. Small coupling: g -C 1 

We can solve for the expansion coefficients in the ma- 
trix ODE ( 56 ) using Standard perturbation theory in the 



small coupling regime. This approximation will be valid 
provided that the mode in question has go -C 1. We first 
expand the matrix A in a power series: 

A = A^+A^g + --- (67) 

and define the following perturbative expansion of c: 

c = c(°)+c< 1 )+-, c<°» = [l0-f 
Then, will satisfy the ODE 
d 



it 



(68) 



dg 



,(D = A (o) c (o) ) c (% ) =0 . 



(69) 



For n m, the elements of the A matrix are given by 



(m — n) J _ 



6 (m — n) 



. (70) 

where the eigenfunctions correspond to the ordinary 
Schrodinger simple harmonic oscillator ( 23 ) . The matrbt 



element is only nonzero if the quantum numbers differ by 
0, 2 or 4, which implies that the only non-zero elements 
of are an d Evaluating the expansion coef- 
ficients in the g — > limit, we obtain the final quantum 
state of the system to be 



I*) 



1 



;sq- 



7.90 



1 



2|2) SQ + ^=|4}sq ) +0(gi), (71) 



9 



which yields the large scale power spectrum 



:1 IA o(— 

2 k* y k^ 



(72) 



where we have made use of the parameterization g$ — 
k/h. 



2. Large coupling: g 2> 1 

It is also relatively straightforward to solve the ma- 
trix ODE (56) in the large coupling regime g ^S> 1. The 



method is similar to the small coupling procedure de- 
scribed above, but now the leading order contribution to 
the coefficient matrix A W is obtained by using the large 
g expansions of the eigenfunctions (51 ) and eigenenergies 
(52) in the matrbc elements (57). Then, assuming the 



same form of the zeroth order solution as before ( 68 1 , 
the first order expansion coefhcients obey 



dc 



(i) 



dg 



5exp [2ia{gl -g 2 )} 
384 5 3 



and 



den 



d) 



dg 



4(n + 1) exp [\ian{n + 2)(gg - g 2 )} 
n 2 (n + 2) 2 (n + 4)(n-2) ff 3 



(73a) 



(73b) 



for n — 4, 6, 8 . . . (with all the coefhcients with n odd 
idcntically equal to zero). Here, we have defined a — 
M+/ 'H. As before, we assume that the mode is in the 
ground state at the beginning of inflation c« (<7o) = 0. 
These ODEs are simple to integrate in terms of Gamma 
functions, but tend to result in long expressions that we 
do not reproduce here. We will compare the results of 
this perturbative analysis to numeric simulations in the 
next subsection. 



D. Numerical solutions 

Note that the perturbative results of the previous sub- 
section can only be used to find final quantum state of 
modes with g -C 1: To obtain the superhorizon be- 
haviour of modes with g$ > 1 we must solve ( 56 ) in the 



transition regime where the coupling is ncither small nor 
large. To do so, we must solve the matrix ODE numeri- 
cally, which involves truncating the infinite dimensional 
system. We will neglect all eigenfunctions with n > n ma x, 
which makes c into an (n max + l)-dimensional vector and 
A into an (n max + l)-dimensional square matrbc. Note 
that A is still anti-Hcrmitian after truncation. 

It is computationally convenient to transform to a new 
time coordinate 



r{9), 



dc 



Bc, B 



dg 

dr 



(74) 



and introduce an evenly spaced r-lattice 
Tj=T -jh, T Q = T(g ). 



(75) 



Here, h is the timestep associated with our numerical 
scheme. Our particular choice for the new time coordi- 
nate is 



T{g) = 2\ng + g 2 



(76) 



For small g, r will be proportional to the cosmological 
proper time t. On the other hand for large coupling we 
will have t m g 2 , which is consistent with the time de- 
pendence of the large g perturbative solutions discussed 
in §IV C 2| 

We write the values c and B at a given lattice point as 



cj = c(s(t,)), B j =B(g(r j )). 



(77) 



A forward-Euler numerical stencil for the solution of ( 74 ) 
is 



/iB,c,- 



0(h 2 



while a backward-Euler stencil is 



c ?+ i = Cj - hB j+ iC j+ i + O (h 



(78) 



(79) 



Taking the average of the forward and backward stencils 
and making use of B J+ i = Bj + O (h) gives 



(I + ±hB j)c j+ i = (I - lhBj) Cj + 0(h 2 ), 



(80) 



where I is the identity matrix. Dropping the error term 
gives our numerical stencil: 



U r 



U, = (I + ihBj)- 1 ^ - i/iBj). (81) 



The advantage of this numerical scheme is that the evo- 
lution operator is automatically unitary U^Uj = I since 
Aj (and hence B j) is anti-Hcrmitian; hence the norm of 
c is preserved: 



,t 

-j+V 



(82) 



This also implies the scheme is unconditionally stable. 
The disadvantages of the scheme are that the global er- 
ror is linear in the stepsize h, and one has to solve the 



linear system (80) for c J+ i at each timestep. For these 
reasons, the method is relatively computationally expen- 
sive to implement. On balance, we find that the unitarity 
and stability of the scheme are worth the additional nu- 
merical overhead. 

A key elemcnt of our numerical analysis involves the 
efficient computation of the a nm matrbc elements (57). 



The Mathieu functions involved in the eigenfunctions ^„ 
are notoriously expensive to calculate numerically, so we 
employ the following strategy: For g < 10~ 2 or g > 
10 2 , we calculate a nm using series expansions [12] of the 
integral 



d 

1 dg 



sm 2 (g 1 / 2 y) 
2ff 



V n dy. 



(83) 
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FIG. 4: Probability |c„| 2 of finding a given Fourier mode in 
the m?" energy eigenstate as a function of the polymer cou- 
pling g. In this example, the mode was assumed to be in 
the vacuum state at an initial time characterized by go — 20; 
i. e. fc/ao-M* = 20. (Note that time runs right to left in this 
plot.) We have selected M* /H = 1. The solid lines are the 
results of a numerical simulation which retains the first nine 
energy eigenstates (n — ... 8) while the dash-dot lines are 
perturbative solutions obtained in Sec. |IV C| 



For 10 -2 < g < 10 2 , we numerically calculate the integral 
at ~ 30 sample points with equal logarithmic spacing 
and then use cubic spline interpolation to deduce the 
integral at other g values. We use a similar combination 
of series expansions and spline interpolation to efficiently 
calculate the energy eigenvalues e n . 

In figure |4j we give an example of the output of our 
numerical code. We note that the numeric and pertur- 
bative solutions closely match for g 3> 1. We also see 
that the expansion coefficients become constant in the 
jC 1 limit. The transition between the two asymptotic 
behaviours occurs for 0.1 < g < 1. 



that for M^/ H > 4, the oscillations appear to be sinu- 
soidal. 

Using perturbation theory, we have already derived the 
large scale limit (72 ) of the polymer power spectrum. We 



can also understand the small scale oscillations of the 
powers spectrum under certain assumptions. Examining 
our numerical results for the evolution of a given mode in 
figure [3J we see the expansion coefficients closely follow 
the perturbative prediction up to some transition epoch, 
and then are roughly constant on large scales. A crude 
approximation to this behaviour is to assume that the 
g — > limit of the expansion coefficients is just given by 
their perturbative values at the transitional epoch g = 
t/tr- More concretely, we can estimate the high go power 
spectrum by evaluating the solutions to ( 73 1 at g = g tr 



and then making used of ( 65 1 and ( 66 1 . 

We can obtain a particularly simple result if we restrict 
our attention to situations where the argument of the 
exponential in (73) is rapidly varying; i. e., M* 3> H. We 
obtain 



2M* 
H 



5y/2 H 
768 M+gi 

fc 2 
fc 2 



K5) 



(84) 



where this expression is only to be applied when fc fc*. 
The form of this approximate power spectrum is con- 
sistent with the simulations results presented in figure 
[5] The deviation of the polymer power spectrum from 
the Standard result is a sinusoidal function of fc 2 on 
small scales whose amplitude is inversely proportional to 
MJH. 



Unfortunately, the expression (84) is of little quanti- 
tative use without knowing the value of the transition 
epoch g t r- However, we can use the functional form to 
motivate a fitting formula for our numerical results that 
is valid on small scales and for M* 3> H. We find that the 
following expression does a reasonable job of reproducing 
simulation results for fc/fc* > 2 and M+/H > 50: 



Vi(k) 



-V 

2tt J 



H 

4M^ 



2M* 
H 



- 1 



(85) 



THE POLYMER POWER SPECTRUM 



VI. OBSERVATIONAL CONSEQUENCES 



The numeric simulations described in §IVD| can be 
used to calculate the power spectrum V^k) using (65) 
and (66). Results for various choices of M+/H are 
shown in figure [5j We see that large scale modes with 
fc -C fc* recover the familiar V$(k) ~ (H /2ir) 2 result from 
Schrodinger quantization, while small scale modes with 
fc 3> fc* exhibit an oscillatory power spectrum. Notice 



In this section, we use the polymer power spectrum 
derived above to calculate cosmic microwave background 
(CMB) angular spectra and the present day (linear) mat- 
ter power spectrum. Our goal is not a detailed compari- 
son to observations, rather we seek to gain a qualitative 
understanding of the polymer effects and an indication 
of whether or not they may be observable. 
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FIG. 5: Power spectrum of polymer perturbations generated during inflation for various values of M*/H. The righthand panel 
illustrates that the power spectrum is an oscillatory function of g$ = k 2 /k 2 for k 3> fc* and that the deviation of the spectrum 
from the Schrodinger result has an amplitude that decreases as M+/H increases. 



A. CMB Angular Power Spectrum 

The CMB angular power spectrum provides the 
highest-quality dataset in modern cosmology, currently 
best-constrained by the WMAP probe [15] ■ The angu- 
lar power spectrum of the temperature auto-correlation 
is given by 



Cti — 



V(k)\A T1 (k)[ 



dk 

T 



(86) 



where Aj<i(k) is the photon temperature transfer func- 
tion, evaluated at the present epoch. Equivalent forms 
hold for the -E-mode polarisation and for the cross- 
correlation. The transfer functions are typically recov- 
ered numerically from a Boltzmann code such as CAMB 
[H] or CLASS [15] . However, we can first gain insight by 
focusing on the large scale, small-Z region of the temper- 
ature auto-correlation, in which 

A T1 (fc) « -~ji[k(r) - 77dcc)] = ~^ji{kx ) (87) 

where rj^ oc is the conformal time at decoupling, r] that 
at the current epoch, and ji is the sphcrical Bcssel func- 
tion of order L This approximation is valid in Einstein-de 
Sitter universes and serves as a reasonable approximation 
for ACDM models with the integrated Sachs- Wolfe effect 
neglected. For a scale-invariant (Harrison-Zel'dovich) 
primordial power spectrum, this transfer function pro- 
duces the Sachs- Wolfe plateau l(l + l)Ci = constant, so it 
could be expected that the polymer quantized primordial 
power spectrum would produce an approximate Sachs- 
Wolfe plateau with oscillations imposed upon it. Sincc 
the polymer power spectrum tends towards a sinusoid for 
high wavenumber we expect the impact on smaller scales, 
and hence higher multipole numbers, to diminish. 



This can be demonstrated explicitly for M+/H > 5 and 
k* <C /c. Across most of the region of integration, k > fc*, 
and the approximate form of the power spectrum ( 84 ) is a 



reasonable approximation even for 1 — 2, and improving 
with increasing l. Consider the fractional change from a 
Harrison-Zel'dovich signal, 



Ci 



C h - C} 



HZ 



H 
4M* J 


sin 


~2KU 
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\kl X ) 




Mkxo)] 2 f 


J 





This equation has an exact solution given by a com- 
plicated combination of hypergeometric and gamma 
functions, inducing the expected oscillating Sachs- Wolfe 
plateau, with the oscillations decaying rapidly. These 
oscillations are of the order of 2% for M+/H = 5. 

The full situation can be solved numerically using the 
transfer functions produced by the CAMB code and ac- 
curate fits for the polymer power spectrum. We take 
fc t = 5x 10~ 4 Mpc -1 and consider the cases M+/H = 1 
and M+/H = 8. We employ a flat WMAP7 concor- 
dance background with h = 0.704, fl b h 2 = 0.02253, 
Q c h 2 = 0.1122, and an amplitude for primordial per- 
turbations of A* — 2.48 x 10 -9 . However, for simplicity 
we take a Harrison-Zel'dovich primordial power spectrum 
with ns = 1; the qualitative results and the amplitude 
of the fractional shifts from the inflationary case are un- 
changed by the spectral tilt. 

Figurerashows the integrands V{k) | A-rz (fc) | 2 /k for the 
inflationary and polymer models. It is clear that this 
choice of k+ provides the best chance of a large impact 
on the CMB, since the peak in the modulation aligns with 
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the peak in the transfer function at low multipoles. 6 

This setup yields the temperature auto-correlation 
Cttu the E'-mode polarization auto-correlation Ceei 
and the cross-correlation Ctei between these, plotted in 
Figure [7] 

The left panel shows the power spectra assuming that 
the polymer power spectrum has the same amplitude A+ 
as that from inflation. However, since in both models 
this amplitude is necessarily a free parameter this would 
be renormalised by the data. The right panel show an 
amplitude modified such that the deviations from the 
inflationary signal vanish for higher L Since WMAP is 
cosmic-variance limited up to the second acoustic peak, 
and Planck is expected to be cosmic-variance limited up 
to l w 2000, this would resemble the best-fit model. 

Figure [8] plots the fractional differences from the infla- 
tionary model, again both assuming the amplitudes to 
be the same (left), and adjusting it to match the higher-Z 
signal (right). This plot shows the form of the deviations 
from the inflationary prediction more clearly. 

The CMB angular power spectra for the model with 
M* /H = 8 are indistinguishable by eye from the in- 
flationary spectra, with a maximum deviation at the 
quadrupolc of ~ 1%, and an asymptotic deviation of 
~ 0.1% even when the amplitude is not renormalised. 

The more extreme model with M+/H = 1 causes os- 
cillations in the Sachs- Wolfe regime, with maxima at of 
the order of ~ 10% at l « 10 in each of the correlations, 
and l w 6 for the cross-correlation. However, these lie 
within cosmic variance around the Harrison-Zel'dovich 
signal and would therefore be extremely difficult to dis- 
tinguish even for this extreme case. As a result it will be 
impossible to recover statistically-meaningful constraints 
from the angular power spectra of the CMB; conversely, 
the lack of such observed oscillations does not imply that 
the polymer quantized models are ruled out. 

It would be interesting to study the primordial non- 
Gaussianity induced by the polymer quantized fluctua- 
tions; if these are of a charaeteristie form then it is en- 
tirely conceivable that there is a greater observable im- 
paet on the CMB bispectrum than on the angular power 
spectrum. We leave this issue to a future study. 



B. Matter Power Spectrum 

The other obvious observable to consider is the matter 
power spectrum 

P M (M) = ^tW)WM)| 2 - (89) 



It also implies that we might alternatively maximize the impaet 
by either choosing fc* to align with a peak in the signal at l = 6, 
or else to align with a multipole corresponding to the top of the 
first acoustic peak. 
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FIG. 6: The integrands of C; when l = 2 for M+/H = 1 
and M+/H = 8. Plotted are the integrands for Ci and for 
ACi/Ci, in comparison with that for the Harrison-Zel'dovich 
spectrum. 

The observable regions of the matter power spectrum 
are in the regime k ^ fc*, implying that for M* /H > 5 
the approximate form of the power spectrum can be used. 
The fractional difference between the matter power spec- 
trum in the polymer quantized model and a Harrison- 
Zel'dovich model is then 

AP M (fc) = P M (fc) ~ fM,Hz(fc) 

P M (k) P M ,Hz(fc) 

H \2M„ f k 2 \1 H , , 
« s n -(ts--1 • 90) 

We then expect sinusoidal oscillations around the 
Harrison-Zel'dovich of order H/4M+, or of the order of 
5% for M+/H 5, regardless of the value of fc*. The 
errors on the observed matter power spectrum are cur- 
rently of the order of 2% [35] , implying that a detection 
may in principle be possible for M+/H < 10. While with 
current observations it will most likely not be possible to 
observe the oscillations imprinted on the matter power 
spectrum, should M+/H be low enough it is entirely pos- 
sible that upeoming missions such as Euclid will observe 
them. 

In Figure [9] we plot the predieted matter power spec- 
trum for the Harrison-Zel'dovich model, and for the poly- 
mer quantized models considered for the CMB. We em- 
ploy an accurate fit for the primordial power spectrum 
rather than the approximation employed above. In con- 
trast to the CMB angular power spectrum, where the 
oscillations induced by the polymer quantisation are in- 
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FIG. 7: The CMB angular power spectra for polymer quantized models with M+/H = 1 (blue) and M*/ H — 8 (black), both 
with fc* = 5 x 10 _4 Mpc _1 . Left: Equal amplitudes for inflationary and polymer models. Right: Polymer model amplitudes 
adjusted for high l. The grey region shows cosmic variance and the model with M* /H — 8 is indistinguishable by eye from the 
inflationary case. 
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FIG. 8: The relative shift in the CMB angular power spectra ACi/Cl = (C; — C; HZ )/C; HZ for polymer quantized models with 
M*/H — 1 (blue) and M+/H = 8 (red), both with k+ = 5 x 10 _4 Mpc _1 . Left: Equal amplitudes for inflationary and polymer 
models. Right: Polymer model amplitudes adjusted for high L The grey region shows cosmic variance. Note that the rapid 
oscillations in the primordial power spectrum tend to cancel themselves out in the angular power spectra for larger M*/H. 



tegrated out except on very large scales, the matter power nonlinear smoothing 



spectrum retains these. Figure 10 focuses on the baryon 



acoustic oscillations, plotting the oscillations around a 
smoothed spectrum. We use a modification of the pro- 
cedure employed by, for instance, [461 147) , and employ a 



Pm 



olh 



(k) 



Pu{k) 



olh 



(k) 



1 



exp 



V2 



-1, (91) 



with Snl = 4.47 to model the damping of the baryon 
acoustic oscillations due to nonlinear processes. Plot- 
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FIG. 9: The matter power spectrum for a Harrison-Zel'dovich 
model (orange), and polymer quantized models with M+/H = 
1 (blue) and M* /H = 8 (black), both with k* = 5 x 
lGr 4 Mpc _1 . 
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FIG. 10: Baryon acoustic oscillatioris for the Harrison- 
Zel'dovich power spectrum and for the polymer power spec- 
trum with M* / B = 8 and K = 5 x 10~ 4 Mpc~ 1 . This case 
should be measurable with the BOSS/CMASS survey. 



ted are the envelopes of the oscillating spectrum rather 
than the spectrum itself, and we neglect the case with 
M* /H — 1; the spread would virtually fiil the plot. Such 
an extreme model is therefore strongly disfavoured by 
present data. Comparison with Figure 18 in [35] suggests 
that while SDSS-II would be unable to detect the case 
with M* /H — 8, it is on the edge of detectability with 
BOSS/CMASS. The upcoming Euclid probe 03] will be 
able to constrain the ration M+/H much more tightly. 



VII. DISCUSSION 

In this paper, we have considered the polymer quan- 
tization of a massless scalar field in an exactly de Sitter 
inflationary cosmology. We have pursued the approach 
of Fourier transforming the field at the classical level, 
which reduces the system to a collection of decouplcd 
simple harmonic oscillators, and then numerically solv- 
ing the resulting time dependent Schrodinger equation 
for the state which closest resembles the Bunch-Davies 
vacuum at the start of inflation. We hence calculate the 
spectrum of primordial perturbations, which recovers the 
scale-invariant Harrison-Zel'dovich result for k less than 
a pivot scale k+ corresponding to a mode with physical 
wavelength 27t/M* at the beginning of inflation [c/. (61 )]. 
Here, M* is the polymer energy scale. For modes with 
k > k*, we find that polymer effects impose oscillations 
onto the Standard result of amplitude oc H /M*. Numer- 
ical plots of the polymer power spectrum are shown in 
figure [5] and semi-analytic fitting formulae are given by 



1 



1 k 

2fc? 
H 

c 

4M t k 



k < fe*, 

2M*/fe 2 \1 fc»fc*, 
H \kl )y M,» H, 

(92) 

where "Phz = (H /2n) 2 is the Harrison-Zel'dovich result. 
These are the principal results of this paper. 

We have also calculated various CMB power spectra 
using the polymer power spectrum. We have seen that 
the largest effect is for low l multipoles, and even for 
H ~ M+ the magnitude of the effects are within the cos- 
mic variance uncertainty. It is hence unlikely that the 
kind of oscillations predicted by this work could be di- 
rectly observed or constrained in CMB angular spectra. 
However, the polymer effects on the present day mat- 
ter spectra are more pronounced, raising the possibility 
that future probes of the baryon acoustic oscillations can 
constrain H/M+ < 0.1. 

We are prevented by performing a more detailed com- 
parison to observation by the fact that our calculations 
have been limited to exactly de Sitter cosmologies. How- 
ever, we expect our results to go through for slow-roll 
inflation: The characteristic inverse timescale associated 
with polymer effects is given by the logarithmic deriva- 
tive of the polymer coupling r~ i ?mer = g/g = -H. On 
the other hand, the characteristic inverse timsescale asso- 
ciated with changes in the background geometry in slow 
roli is the logarithmic derivative of the Hubble parameter 



H /H. Hence we find 



^polymer 
^slow-roll 



H 



= e < 1, 



(93) 



where e is one of the Standard slow-roll parameters. 
Hence, we nai'vely expect that our polymer power spec- 
trum results to be valid in the slow roli approximation 



with the H appearing in ( 92 ) being interpreted as the 
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Hubble factor at horizon exit for a given mode, but more 
detailed calculations are warranted. 

It is interesting to revisit the model presented in Ref. 
[35] . which found that the polymer quantization of a 
homogeneous scalar could drive quasi-dc Sitter inflation 
with Hubble parameter 



H 2 = 



For this model, we then have 



H 

W* 



1 E, 



inf 



\/6 -Mi 



1.6 x 10" 



pi 



E; 



inf 



10 16 GeV 



(94) 



(95) 



We see that for reasonable choices of the inflationary en- 
ergy scale, H/M+ will be small for these models. This 
implies that polymer effects on the perturbation spec- 
trum will be rather small and likely unobservable. Fur- 
thermore, the slow roli parameters for this model are 
exponentially small [35] . which implies that generated 
nuctuations will be very close to the Harrison Zel'dovich 
result. WMAP data appears to disfavour a purely scale- 
invariant spectrum [35] , which suggests some tension be- 
tween polymer-driven inflation and current observations. 

Finally, we remark that our technique of Fourier trans- 
forming a scalar field at the classical level and then quan- 
tizing seems to provide a novel avenue for including high 
energy/small distance modifications into inflationary cos- 
mology. It is reasonable to assume that different effects 
(such as modified uncertainty relations) would lead to 
different classes of time-dependent potential appearing in 
the Schrodinger equation ( 40 ) , and hence different power 



spectra. We will report on such models in the future. 
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Appendix A: Polymer quantum mechanics 

In this appendix, we describe the polymer quantiza- 
tion of a particle moving in the x direction under the 
influence of a potential V(x). In {IV we will apply the 



same techniques to quantize the amplitude of a single 
Fourier mode of a massless scalar during purely de Sitter 
inflation. 

The Hamiltonian of this simple system is 



1 

2ro 



P 2 + V(x), 



(Al) 



where p is the momentum canonically conjugate to x, and 
we have the Poisson braeket {x, p} = 1. In polymer quan- 
tum mechanics, it is assumed that the quantum state of 



the system can be expressed as 



h/>> = E 



3=-oo 



Here, \xj) are eigenstates of the position operator 



(A2) 



(A3) 



and the c/s are expansion coefficients. Note that we 
will assume that all the time dependence of the sys- 
tem is carried in the state veetor = \ip)(t) (i. e., we 
take the Schrodinger picture), which means the expan- 
sion coefficients will in general depend on time Cj = Cj (t) . 
The spectrum of the position operator {xj } consists of a 
countable selection of points from the real line K which 
is analogous to the graphs covering 3-manifolds in LQG. 
A critical difference between polymer and Schrodinger 
quantizations is that in the former, the position eigen- 
states are normalizable: 



(xi\xj) = 5i 



(A4) 



A translation operator can be defined by it's aetion on 
position eigenstates: 



U x \xj 



A); 



(A5) 



that is, U\ converts a position eigenstate with eigenvalue 
Xj into an eigenstate with eigenvalue Xj — X. One can 
easily confirm the commutator identity 

[x,U x ] = -XU x , (A6) 



by aeting the lefthand side on an arbitrary state ( A2 ) . 

In Schrodinger quantum mechanics, the relationship 
between the translation and momentum operators is 
i. e., p is the generator of infmitesimal trans- 



Ih 



lations. Another way to state this is: 



.au x 



0X 



(A7) 



A=0 



Hence in order to define a momentum operator, U\ must 
be differentiable at A = 0. It is fairly easy to see that this 
is not the case: for example, consider the expectation 
value of U\ in a position eigenstate: 



(xi\U\\xi) 




(A8) 



In is elear that the matrbc representation of U\ in the 
basis (and hence U\ itself) fails to be differentiable 
at A = 0. Therefore, we cannot define a momentum 
operator in the polymer setup. 

However, we can define an effective momentum opera- 
tor by creating a finite difference stencil of the definition 
(|A7|: 



PA» 



2A* 



U\< - U 
2iA* 



(A9) 
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(This is not the only way to define p\ t , see footnote [4] 
for more details.) The fixed width of our finite difference 
stencil A* plays a crucial role in this quantization schcmc: 
As mentioned above, it defines and energy scale M* above 
which Schrodinger and polymer predictions will diverge. 7 
Practical calculations in polymer quantum mechanics in- 
volve mapping p i— >• p\jn the operator version of the 
classical Hamiltonian (A10): 



1 

2^n 



PK 



fV(i) 
1 

8mXl 



(2 - &2X t 



Ul) + V(x). (A10) 



Armed with this operator representation of the Hamilto- 
nian, we can attempt to solve the Schrodinger equation 



(Ali) 



and hence determine the time evolution of observable fea- 
ture of our system. 

Now, the graph {xj} involved in our decomposition of a 
generic state vector \ip) is meant to be an arbitrary selec- 
tions of point on the real line, not necessarily a regularly 
spaced lattice. However, we find that the Hamiltonian 
tends to "super-select" regular lattices in the following 
sense: We say that a state vector is in the subspace 7^°^ 
of the entire polymer Hilbert space H po i y if it can be ex- 
pressed as 



An arbitrary state in this basis is expressed as ip(p) 
(p|V>)- It follows that 



oo 



c j\ x j) 



w, p) 



oo 



/ j -J t~ J I T\-ir/ / , ^J' 

j — — oo j — — oo 

(A15) 

In Schrodinger quantum mechanics, \p) would be an un- 
normalizable momentum eigenstate with the above sums 
replaced by integrals; however, the lack of a momentum 
operator in polymer quantization prevents us from mak- 
ing a similar interpretation here. Nevertheless, we call 
ip(t,p) the momentum-space wavefunction. It is easy to 
deduce the aetion of fundamental operators in this rep- 
resentation: 



W,p). (A16) 



In this basis, the Schrodinger equation |All| becomes a 
diffcrcntial equation: 



2m\l \ dp 



i>(t,p). (A17) 



Whether or not this is easier to solve than the difference 



equation implicitly defined by (Ali ) depends on the form 
of the potential. 

Let us now restrict ourselves to a particular super- 
selected subspace Hp° ly . An identity operator on such 
a subspace is 



= <T c j \x + 2X i< j) & eH x p ° oly , (A12) 



where A* is the fixed scale we used in the definition of 



Px ii (|A9|). It then follows that if \ip) is in "Hp° ly then the 



Hamiltonian (AlOl aeting on 



is also in ^° ly ; 



fc rt poly 



(A13) 



Hence, the Schrodinger equation (Ali I and the as- 



sociated eigenvalue problem can be solved in each of 
the super-selected seetors separately. In more practical 
terms, this means that it sufhcient to assume state vec- 
tors of the form (A12) in ( All[ ), which translates the 
problem into an infinite dimension matrix ODE for the 
expansion coefficients c = [c,- (t)]. 

The solution of such a problem can be unwieldy, so we 
introduce an alternative representation. Let us define a 
new basis via 



\P) = 



OO 

E 



'>,•>, (p\ XJ ) = 



(A14) 



7 The precise relationship between A* and M* depends on the 
dimensions of x: If x is a length, then M* = l/A*. 



7T/2A* 



-tt/2A« 

It is easy to confirm that if | 
then 



dp\p)(p\. 



e Klx y and 



(A18) 



(Z ']-/ x 

c rL poly' 



r/2A, 



(</#> = (0iiiv^> = 



r(t, P Mt,p)dp, (A19) 



-7T/2A* 



which serves to define the inner produet on %poi y in the 
momentum representation. 

Finally, we note that if we restrict ourselves to states 
in "Hp° ly , we can take p e [— 7r/2A*,7r/2A*]. From the 



definition ( |A15 ) and xj = xq + 2A*j, the momentum 
space wavefunction must satisfy the following boundary 
conditions: 



il>(t, -7T/2A*) = e^°/ A ^(i, 7T/2A,). 



(A20) 



A specific choice of lattice offset will result in a specific 
boundary condition: For example, if we choose Xq = 0, 
the wavefunction will be periodic with period 7r/A*, while 
the choice xq = A* results in an anti-periodic wavefunc- 
tion. This approach to selecting boundary conditions will 
have physical ramifications; in the case of the polymer 
simple harmonic oscillator, the energy eigenvalues will 
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depend on Xq [501 IH]- I n this paper, we take a different 
approach by demanding ip(t, — 7r/2A+) = tp(t, n/2\±) = 0. 
In, §IV A| we argue that such a choice is required to en- 
sure unitary evolution with a time-dependent A. Here, 
we see a separate justification: This is the only choice of 
boundary condition that will ensure physical observables 
do not depend on x ; i. e., the choice of regular lattice in 
which we perform calculations. 



the start of inflation. Then, 



7 <^cnd ^start * T 

K* = M*. 

^now ^cnd 



The number of e-folds of innation is 



N = ln 



&cnd 
^start 



(Cl) 



(C2) 



Appendix B: Scaling formulae 

In an FRW universe, all physical or observable quan- 
tities must be invariant under dilations of the spatial co- 
ordinates, or equivalently, re-scalings of the scale factor 
a. It is useful to make note of how the various quantities 
defined above transform under a uniform scaling of the 
spatial coordinates by a constant factor t: 



o. 



6 k ->■ ^ 3/2 0k, 



7T -j» £ 7T, 
r 3/2 



k £ _1 k, 



7T k . 



(BI) 



We use these relations in S IV to define field amplitudo 
translation and effective polymer momentum operators 
that transform correctly under dilations x — > £x. 



Appendix C: Derivation of fc* 



The radiation-like component of the universe's total den- 
sity evolves adiabatically after inflation. Equating the 
entropy in the radiation fluid now with the entropy just 
after the end of inflation yields 



Ganda. 



3 rpS 

end end 



'now»now-'iiowi 



(C3) 



where Q is the effective number of relativistic species in 
the universe at a given epoch, and T is the temperature. 
We can use the Stepan-Boltzman law to relate the tem- 
perature to the total density of radiation at the end of 
inflation: 

Pend = ^^endTT^bJcLj, (C4) 

Where k\, is Boltzmann's constant. Just after reheating, 
the universe is radiation dominated so p cno \ is actually 
the density appearing in the Friedman equation: 



Tj2 Tj2 



Pcnd 

3Mf l 



E inf 

3M|, : 



(C5) 



Most of our calculations and results are defined in 
terms of a quantity fc* which is the present wavenum- 
ber of a Fourier mode whose physical wavenumber was 
equal to M+ at the start of inflation. In this appendix, 
we derive an explicit formula for fc* in terms of familiar 
inflationary parameters. Define a now as the current value 
of the scale factor, a cn( j as the value of the scale factor at 
the end of inflation, and a sta rt to be the scale factor at 



where H is the Hubble parameter during inflation and 
-Eirif = Pend ^ s inflationary energy scale. Assembling 
these results, we obtain 



ir 1/2 k h T now Q 



1/3 
now , 



3 • 30V4 



^end 



(C6) 



Using T now = 3.94 K and G now = 3.04, we obtain (61 1 
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